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At arbitrary temperature T, we solve for the dynamics of 
single molecule magnets composed of three classical Heisen- 
berg spins either on a chain with two equal exchange con- 
stants Ji, or on an isosceles triangle with a third, different 
exchange constant J2. As T — > 00, the Fourier transforms 
and long-time asymptotic behaviors of the two-spin time cor- 
relation functions are evaluated exactly. The lack of trans- 
lational symmetry on a chain or an isosceles triangle yields 
time correlation functions that differ strikingly from those on 
an equilateral triangle with Ji = J2. At low T, the Fourier 
transforms of the two autocorrelation functions with Ji 7^ J2 
show one and four modes, respectively. For a semi-infinite 
J1I J\ range, one mode is a central peak. At the origin of this 
range, this mode has a novel scaling form. 



I. INTRODUCTION 

Recently, there has been a substantial interest in the 
physics of magnetic molecules, or single molecule mag- 
nets (SMM's). Such studies are important 
both for basic scientific reasons and for possible techno- 
logical applications. Among the smaller SMM's are 
dimers of {S = 1/2) and of Fe3+ (5 = 5/2V [|jlO| 
a nearly equilateral array of three spins, and 
an isosceles triangle of Gd^+ {S = 7/2) spins (Gd3). 
Non-metallic variations of OL-BaRuOa with three Ru^ 
{S = 1) ions (Ru3) might behave as three-spin chain 
SMM's. 

Each SMM consists of a small number of paramag- 
netic ions surrounded by non-magnetic chemical ligand 
groups, and is large enough that the magnetic interac- 
tions between the ions in different SMM's within a crys- 
tal are negligible. Hence, measurements performed on 
macroscopic samples just probe the magnetic interactions 
within the individual SMM's. Large single crystals with 
long-range SMM packing order are suitable for inelastic 
neutron scattering experiments. Measurements at the 
appropriate deviations from the Bragg wave vectors for 
the particular SMM crystal structure probe statistical 
ensembles of the Fourier transforms of the two-spin time 
correlation functions involving the spins in each SMM. 

The magnetic interactions between the ions in these 
SMM's can often be described by the isotropic Heisen- 
berg model. One expects that for ions such as Gd^+, 
with S — 7/2, the classical version of the Heisenberg 
model captures the essential features of the dynamics at 
not too low temperatures T. For a dimer, comparisons 
of the classical and S — 1/2,5/2 quantum behaviors of 



the dynamics supported this expectation. fl^jQ For the 
equilateral triangle, such comparisons were only made 
for T 00. [T^ In addition, exact and numerical results 
were presented for the classical dimer |l^Jl^,^ , the four- 
spin ring pT[ | , and the A^-spin equivalent neighbor model, 
which includes the equilateral triangle, 

For three spins on an equilateral triangle, with equal 
Heisenberg exchange constants Ji , the spin sites on each 
triangle are translationally invariant. Spins on an isosce- 
les triangle have one additional exchange constant J2 ^ 
Ji, and this simple translational invariance is absent. 
Dynamical measurements on Gd3 and modified composi- 
tions of Ru3 could help to uncover the effects of this lack 
of translational invariance inside each SMM, and might 
aid in our understanding of more complicated systems 
with multiple magnetic interactions. 

At low T, the A^ spins in the equivalent neighbor model 
oscillate in a single mode. [|2^ Here we show that the ab- 
sence of translational symmetry inside each isosceles tri- 
angle generally introduces three additional low-T modes 
at tunable frequencies depending upon 7 = J2/Ji- For 
a semi-infinite range of 7 values, one of these additional 
modes is a central peak. In addition, the dynamics of 
the spins on the endpoints of the three-spin chain with 
J2 = are qualitatively different from the dynamics of 
the spin at the chain center. The mode tuning parameter 
7 makes the dynamical behavior of the spins on an isosce- 
les triangle remarkably different from that present in any 
SMM system studied previously, including the four-spin 
ring. 

The structure of the paper is as follows: In Sec. || 
we define the model, calculate its partition function, and 
present the exact time evolution of the spin vectors. In 
Sees. [H and IV, we discuss the time dependence and 
Fourier transform, respectively, of each spin-spin corre- 
lation functions. We discuss our results in Sec. ^ 



II. THE MODEL AND THE SPIN DYNAMICS 

We study the Hamiltonian describing three spins of 
unit magnitude, 5'i = |Si| = 1, on a triangle with two 
classical Heisenberg exchange couplings. 



H — — Ji(Si •82 + 82- 83) — J2S1 • S3. 



(1) 



The cases Ji = J2 and Ji 7^ J2 ^ describe equilateral 
and isosceles triangles, respectively. The case J2 = de- 
scribes a three-spin chain with free boundary conditions. 
Equation can be realized if the ring contains either 
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two different lattice constants, as in Gd3, or ions with 
two different spin values. 

We rewrite Eq. (|l|) in terms of the total spin 8 = 82 + 
S13, with S13 = Si + S3, and obtain, up to a constant. 



2 2 

The partition function Z = / (HLi dn,/ATT)e~'^" is 

Jo A^-M 
= — / dxe°''^''^ sinh(2aa;). 



(2) 



(3) 



where /3 = {kBT)~^, dfti is the solid angle element for 
the i-th spin, a = f3Ji/2, 7 = J2/ Ji, x = S13, and s — S. 

To calculate the time-dependent correlation functions 
(Si(i)-Sj(O)), we first solve the classical Heisenberg equa- 
tions of motion for the quantities S2(i) and Si3(i), 

82,13 = ^18243 X S, (4) 

leading to 

S2,i3(i) = C2,i3S -t- ^2,13 [cos(Jis<)x - sin(Jist)y] , (5) 

where s||S, s = x x y, A2 = -A13, C2 + C13 = s, C13 = 
(s2 + ^2 _ i)/(2s), and Aj^ = x^ - Cf^. 

The time dependence of Si (or S3) is obtained from 



Si = (J2- Ji)Si X S13 + JiSi X S, 
leading to 

Sis{t) = '5'iso + AS'iso cos[(Ji - J2)xt - (/)o] 
bi±\tj = — -p:, exp(=F«sJit) 



(6) 



(7) 



+ 



Ci3 

AisASisQ 



2(Ci3 + 
AiaASiso 
'2{Ci3-x) 



exp {T'iiiJis + [Ji — J2)x]t ± i0o} 
exp {=Fi[( Jis - ( Ji - J2)x]t T iM , (8) 



where Si± = Six±iSiy, Siso = C13/2, ASuo = ^13(1 - 
x'^/4:)^/^/x, and (/)o is an arbitrary angle describing the 
initial spin configuration. pi[ 



III. TIME CORRELATION FUNCTIONS 

We study the two-spin time correlation functions 

C,,(t)==(S,(t).S,(0)), (9) 

where (...) = J^^ dq^„ dx J^+_\^sdse~^" ... /i27rZ). 
Since Cu — C33 and C12 = C23, only four Cij are inde- 
pendent. These are constrained by the sum rule 



= C22(t) + 2Cii(t) + 4Ci2(t) + 2Ci3(0- (10) 

After averaging over (/)o , the exact Cij {t) satisfy 

Cii(<) =Io + hit) + hit) + hit), (11) 

Cuit) =h + hit) - hit) - hit), (12) 

C22it) = 4 [/o + hit)] + 2(C2s) - (13) 

Ci2(<) = -2 [h + hit)] + \ [is') - iC2s)] , (14) 



where by setting t* = Jit, we have 

h = lic!,), 

hit)^\iAl,cosist*)), 
l.A^ 



(15) 
(16) 



^3(i) = (^^^((Ci3 -x)2cos{[,s+ (1 -7)x]r} 



+ (Ci3 + a;)' cos{[s - (1 - 7)a;]r }) ). 



(18) 



As T — > 00, the remaining double integrals in Eqs. 
(|l5|)-(|l8|) can be reduced to single integrals, as shown in 
the Appendix. Note that hit) is a constant for 7 = 1, but 
not for 7 7^ 1, leading to dramatic differences between the 
dynamics of spins on equilateral and isosceles triangles, 
respectively. As t ^ c«, the time-dependent i nteg rands 
oscillate wildly and their contributions vanish. 1 23 1 From 
Eqs. ( p^ and (^6|), lim^.t^oo C22 for 7 ^ 1 is identical to 
its value for 7 = 1. Only at finite T docs J2 influence 
C22it). Hence, for T > 1 and T 00, 



^ 1 sin(r) -hsin(3r) 
hm C22(t) -- + 25s (^,)3 



T—'Xi 



.3 = |jln3- 1. 0.07359. 



(19) 
(20) 



From Eqs. (|l|)-(|l|), it then follows that 



T n 1 . sin(r) -Hsin(3r) 

l-Ci2(i)^--^3 + ^4(^:,y^. (21) 



t*>l 



The surprise comes from the behaviors of Cu and C13. 
The long-time approaches to the various asymptotic lim- 
its are characterized by three different powers, corre- 
sponding to the equilateral, the isosceles, and the chain 
cases, respectively. For 7 = 1, Cj^^it) = C22(t), which 
approach their mutual limit as a.s T —> 00. 

For 7 = 0, Cii(t) is dominated by hit), yielding 



lini Cjr'it) 



1 



2 



1 

2t* 



■ sin(t* 



(22) 



which approaches its different asymptotic limit dramati- 
cally slower. 
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For 7 7^ 0, 1, hit) and Isit) decay as t~'^ and oscillate 
with different frequencies. The evaluation oflsit) is quite 
involved. By calculating its Fourier transform and then 
inverting it through integration by parts, we find, 



lini 

T|t*,|l-7|t* 



S3 

2 



{ai + 02 cos[2(l - 7)i*] + as cos(t*)} 



it*) 



*^2 



where (see Appendix) 
1 

0.1 = 



0.2 
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6(1-7)2^ 
1 / 5 
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Note that one cannot take either of the limits 7^1 
or 7 ^ directly in Eq. ( p3| ) , since the expansion is valid 
only when both |7|i* » 1 and |1 -7|i* > 1 are satisfied. 

In Fig. |l], we plot the Cij{t) as T 00 for 7 = 0. Since 
C22{t) and Ci2(t) as T ^ 00 are independent of 7, these 
curves are respectively identical to those for Cii{t) and 
Ci2it) obtained for T ^ 00 in the equilateral triangle. 
These functions each approach their asymptotic limits as 
t~^. On the other hand, for 7 = 0, Cii(t) and Ci3{t) 
oscillate about each other for a long time, approaching 
their mutual asymptotic limit as t^^. 
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FIG. 1. Plots of Cii(t) (solid), Ci2{t) (dashed), Ci3(t) 
(dense-dotted, Ci3(0) = 0), C22(t) (sparse-dotted, C22(0) = 1) 
versus t* = \ Ji\t as T ^ cx) (a = 0) for the chain (7 = 0). 

In Fig. H we compare Cii(t) and Ci3{t) as T ^ 00 for 
7 — ±0.3 and 1.5. For 7 = ±0.3, the leading term in the 
approach to the asymptotic limit arises from 03 in Eq. 
(p3|). For 7 = 1.5, 2|1 — 7I = 1, so there is asymptotically 
only one frequency, equal to that for the chain shown in 
Fig. 1, but the asymptotic limit is approached faster. 



At finite T, the physics of the model is influenced not 
only by 7, but also by the sign of Ji. We henceforth refer 
to the Ji > and Ji < cases as ferromagnetic (FM) 
and antifcrromagnetic (AFM), respectively. 
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FIG. 2. Plots of Cii(f), Cia{t) versus t* = \ Ji\t for T 
(a=0), 7 = 0.3 (solid), -0.3 (dashed), and 1.5 (dotted). 




FIG. 3. Plots of limt-,00 Cij{t) for the chain (7=0) versus 
|q| = /3| Jil/2, for a > (FM) and a < (AFM). 

We obtain the Cij{t) at finite T by direct numeri- 
cal evaluation of the double integrals in Eqs. (pl|)-(p^. 
As an example, in Fig. ^, we plot the linit^co Cij{t) 
for 7 = as functions of |a| = (3\Ji\/2. The FM 
and AFM cases are distinguished by arrows. We recall 
that limt^oo Ci3(i) — linit^oo Cii(t) for all 7 7^ 1. The 
spins are intrinsically unfrustrated. At low T, the FM 
limt^ooCy (i) ~ 1, while the AFM limt^oo Ci2(t) -1. 
In the AFM case, limt_ 
- 0.29 at a 0.9. 



, Cii(f) has a minimum value of 
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IV. FOURIER TRANSFORMS 

The dimensionless Fourier transform (FT) of each de- 
viation dCij{t) EE Cij{t) - Muit^ oo Ci J (t), HI] is 

I T I P+OO 

SC^j{lj) = ^-^ / dte*"*(5C,j(0. (27) 

J-oo 

Since causahty requires SCij{uj) — SCij{~uj), we consider 
only positive values of lo j J\. The Fourier integral is in 
principle elementary, since the time dependence is con- 
tained in simple trigonometric factors, yielding combi- 
nations of 5-functions that allow us to perform one of 
the integrations. We thus obtain single integral repre- 
sentations for the FT's at any T, which can be evalu- 
ated with high precision. Since each integration domain 
is restricted to a finite region in the (s,x) plane, the b- 
functions only contribute to the integrals if w falls within 
specific ranges. Determining the allowed ranges of lo val- 
ues is tedious. Some of the details of this calculation are 
presented in the Appendix. 

A. Infinite temperature 



discontinuity at Ji| = 2, the chain i5Cii(a') is discon- 
tinuous at a;/|Ji| — 1, as indicated by the symbol ~. 
This discontinuity is responsible for the 1/t long-time, 
T — > oo decay of 5Cii[t). 

The isosceles triangle cases are also very interesting, 
as they do not just interpolate between the equilateral 
triangle and the chain. For 7 = ±0.2, the discontinuities 
in the slope of SCuiu) at w/|Ji| = 1,1.4, and 1.6 are 
also discernible. As 7 ^ 0, the slope discontinuities at 
'^/I'^il = 1 approach ±00 as expressed in Eq. (|2^). In 
addition, there is a large discontinuity in the slope of the 
curve for 7 = 0.9 at w/Ji — 0.2, and a resulting large 
peak at w/jJij k, 0.088, arising from /2(w). This slope 
discontinuity is responsible for the first two terms oc 
in Eq. (|3|). 

B. Low temperature 

The 5Cij{uj) at finite T were obtained by numerical 
evaluation of the single integrals defined in the Appendix. 
As T — > cxD, we compared our numerical and analytic 
results. At finite T we also checked the initial value sum 
rules by integrating the FT's in Eq. (p7|). 



T ^ 00 




C0/IJ1I 



FIG. 4. Plots of 5Cii{u) as functions of lij/IJil as T ^ 00, 
for 7 = 1,0 (solid), 0.2 (dashed), -0.2 (dotted), and 0.9 
(dash-dotted). 

In Fig. ^ we plot 5Cii{uj) asT ^ 00 for some I7I < 1. 
The solid curves for 7 = 1,0 correspond to the equilat- 
eral triangle and the chain. For the equilateral trian- 
gle, 5Cii{u}) vanishes as w — > 0,3, and exhibits a single 
smooth peak at w/|Ji| « 1.4385. Although difficult to 
discern in this figure, the curve contains discontinuous 
second derivatives at a;/|Ji| — 1,3 responsible for the 
t-3 long-time, T ^ 00 decay of 5Cii{t). ||2|,|2| 

As r ^ 00, although 5C22{'^) for the chain coincides 
with the equilateral triangle 5Cii{lo) curve, the chain 
5Cii{uj) is dramatically different. In addition to the slope 



1. Equilateral triangle 
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FIG. 5. Plots of \a\-^^'^5Cii{u) versus \a\^'^uj /\Ji \ for the 
AFM equilateral triangle at low T. 

We first discuss the equilateral triangle, 7 = 1. For the 
AFM case. As T decreases, the single peak in (5Cii(ijj) 
grows in amplitude and shifts to lower lj/\Ji\. As shown 
for the equivalent neighbor model, j2^, this behavior 
can be quantified by plotting \a\^'^/'^5Cii{uj) versus lj — 
\a\^/'^LL}/\Ji\. In Fig. ||, the \a\~^/^dCii{uj) curves ap- 
proach the uniform AFM equivalent neighbor model form 
j^o)^ exp(— w^), shown as the solid curve, as T — > 0. 
p2| This frustrated behavior results in a scaling of the 
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time, as C{t) approaches a uniform function of {JT^^^t 
as T ^ 0. 

For the FM equilateral triangle, as T decreases, the 
peak in 6Cii{lu) shifts asymptotically to 3Ji, as all of 
the spins oscillate together. In the equivalent neighbor 
model, as T — * 0, the three-spin 7 = 1 FM curves 
approach a uniform function of {uj/Ji — ^^^3)0, where 
= 3 - l/(3|a|), as shown in Fig. |. H That is, as 
T ^ 0, the peak amplitude tends to a constant value, 
but its maximum position approaches 3Ji linearly in T . 
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FIG. 6. Plots oi5Cii{u)/A3 versus (tj/Ji - ajJ)3a/2^/^ for 
the FM equilateral triangle at low T, where A3 = 8/(3e'^) and 
c^3* =3-l/(3|a|). 



2. Three-spin chain 
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FIG. 7. Plots of 5Cii(w) versus (cj/|Ji! - l)2|a|/5 and 
SC22{uj) versus (a;/|Ji| - l)2|a|/5 -1-0.2, for the AFM chain 
at the low-T values -a = 5,10,20,40,80. Cii{uj) is discon- 
tinuous at (jj/|Ji| = 1, indicated by the ~. 

Strikingly different behavior is obtained for the low-T 
AFM chain. In this case, instead of the single AFM mode 
present in the equilateral triangle, there are four modes. 



the frequencies fii of which approach fli/lJi] = 1,2,3 as 
T ^ 0, as the frequencies of the strongest modes, fii, 1^2, 
become degenerate. In Fig. |^, we plotted dCii{uj) and 
(5C22(w) for these two modes, versus (w/l Ji| — 1)2|q!|/5, at 
various low-T a values. Curves for 6C22{^^) are shifted to 
the right by 0.2 for clarity. We note that 5Cii{ui) for the 
^2 mode has a shallow maximum at a frequency which 
approaches | Ji| from above as T ^ 0. For the fii mode, 
SCii has a large discontinuity precisely at Ji | = 1 for 
all T. By contrast, the 6C22{t^) curves only exhibit the 
fl2 mode, and are smooth and rather symmetric about 
their maxima. Both SCii{u!) for these modes approach a 
uniform scaling function of |a|((jj/| Ji| — 1) as T — -> 0. We 
remark that the low-T scaling behaviors of the dominant 
modes in the AFM chain are similar to that of the FM 
equilateral triangle, since neither of these mode energies 
approaches as T — > 0. 

In addition, we found that the weaker fl^ and 
modes, for which SCii{uj) is peaked at « 3,2, 

respectively, both approach uniform low-T scaling func- 
tions. For the ^3 mode, \a\SCii{uj) approaches a uni- 
form scaling function of |a|(w/|Ji| — 3), and for the 
mode, a^6Cii{u!) approaches a uniform scaling function 
of |a|(w/|Ji| -2) as T ^ 0. Hence, as T ^ 0, the ampli- 
tudes of these modes vanish as T and T^, respectively. 
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FIG. 8. Plots of n^5Cii(Lj) and n^(5C22(t^) versus 
(a;/ 1 Ji I — Lj* ) I a I /n for the VL\ and f22 modes of the FM chain at 
the low-T values a = 5, 10, 20, where oj^ = n— (n— 1) / (2njQj). 
Cii(tj) for the mode is discontinuous at \ = 1, as in- 

dicated by the ~. 

For the FM chain at low T, there are four non- 
degenerate modes for which (5Cii(tj) is peaked at 
ilj/|Ji| K. 1,2,3,5. At low T, the two largest intensity 
modes ^2 and lli, which are peaked at w/Ji w 1,3, re- 
spectively, are pictured for a = 5, 10, 20 in Fig. |[ At low 
T, ^Cii for the ^2 mode is very large for w/l Ji| < 1, but 
drops discontinuously to zero at w/| Ji| — 1, as shown in 
Fig. ||. For the next largest intensity mode, f^i, both 
(5Cii(w) and the larger 5C22(w) exhibit peak positions 
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that approach uj/\J\\ = 3 as T 0, as shown in Fig. 
^. We note that we have multiphed the intensities of 
these curves by a factor of 9, and squeezed the scahng 
variable by a factor of 3, relative to those of SCu for the 
U,2 mode, in order to fit all three sets of curves in the 
same figure. All of the curves shown in Fig. ^ scale as for 
the FM equilateral triangle. As for the AFM chain, at 
low r, a5Cii{ijj) for the 1^4 mode approaches a uniform 
function of (cj/| Ji| — 2)|q!|, and a^5Cii{uj) for the weakest 
mode fis approaches a uniform fimction of Ji | — 5)|a|. 



3. General isosceles triangle 

We now consider the more general isosceles triangle 
cases. We first discuss the AFM cases, Ji < 0. As for 
the chain, for 7 ^ 1, as T is lowered, i5Cii(a') generally 
develops into three or four peaks, which become progres- 
sively sharper. When only three modes are present, one 
of them is a central peak at a; — 0. For non-vanishing 
frequencies, these modes are magnons. As for the three- 
spin chain, their relative intensities at finite T are very 
different, as the two largest are typically a few orders of 
magnitude larger than the third, which is a few orders of 
magnitude larger than the fourth. 

To illustrate the types of low-T AFM behavior, in Fig. 
^, we plot logio[^Cii(a;)] and logio[^C22(w)] versus i^/\Ji\ 
at the AFMJow-T value a = -80 for 7 = -0.1, 2.5. For 
7 — —0.1, 6Ciiluj) exhibits four sharp peaks at w/| Ji| « 
1, 1.2, 2.2, and 3.2. For 7 = 2.5, there are three broad 
peaks in 5Cii{u)), one of which is a central peak at w = 0, 
and the others are centered at ~ 0.6 and 1.2. In 

each case, 5C22{'^) is a single peak at one of the larger 
non- vanishing 5Cii{lo) peak positions. 
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FIG. 9. Plots of the AFM \og-^Q[5Cijiuj)] versus tj/|Ji|, at 
low T {a = -80) for 7 = -0.1, 2.5. The solid (dashed) curves 
correspond to SCuiuj) for 7 = —0.1 (2.5). The dash-dotted 
(dotted) curves are the corresponding 5C22(<^) curves. 

After a careful analysis of many lower-T results, we 
established simple formulas relating the frequencies of 



the modes to 7. 5Cii{u!) contains four low-T AFM mode 
frequencies, ili(7), depicted in Fig. ^ which satisfy 

l]i(7)/|Ji| - 1, for 7<l/2, 

= II-I/7I, for 7> 1/2; (28) 
f^2(7)/|Ji| = 1-27, for 7<l/2, 

= 0, for 7 > 1/2; (29) 
f^3(7)/l>/i| =3-27, for 7<l/2, 

= 2|1-1/7|, for 7>l/2; (30) 
!^4(7)/l'/i| =2-27, for 7<l/2, 

= |1 - I/7I , for 7 > 1/2. (31) 

5C22{^) contains only the low-T FM mode with fre- 
quency ill (7). We note that for 7 > 1/2, il4 and Vli 
are degenerate. We analytically confirmed Eq. (^8|) by 
performing an asymptotic evaluation of the integral rep- 
resentation of 5C22{^)i as shown in the Appendix. 
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FIG. 10. Plots of the low-T AFM mode frequencies 
S7i(7)/jJij versus 7. The triangles indicate low-T scaling. 
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FIG. 11. Low-T plots of 5Cii{u) vs. the scaled frequency 
'|a|^''V|J'i| displaying the AFM f72(l/2) mode. 

For the AFM case with 7 < 1/2, there are four modes 
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at finite frequencies. For 1 7^ 7 > 1/2, fii and ^4 are 
degenerate, and the mode with frequency ^2(1) = is 
a central peak. The two triangles in Fig. O indicate 
special points. For 7 = 1, the four modes are all de- 
generate, SCii{uj) — i5C22(w), and dCii{uj) approaches 
the AFM scaling form shown in Fig. For the low-T 
AFM case 7 = 1/2, 6Cii{uj) exhibits the mode frequency 
fi2(l/2) I J2|/|a|i/2 as T ^ 0, and the shape of this 
mode approaches a scaling function of |q:|^/'^(jj/| Ji|, as 
shown in Fig. |ll|. These SCii{u!) curves scale without a 
corresponding rescaling of SCu^uj) itself. 
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(a/|Ji| 

FIG. 12. Plots at Q = 20 of the FM logi(,[(5Cii(a;)] for 
7 = 1.5 (solid) and -1 (dashed) and log]^Q[5C22(<i')] for 7 = 1.5 
(dot-dashed) and -1 (dotted) versus a;/|Ji|. 

We now consider the FM case at low T. To illus- 
trate the types of behavior found, in Fig. 03, we plot 
logjQ[(5Cii((jj)] and logjQ[(5C22(w)] versus cj/j Jifat a = 20 
for 7 = 1.5, —1. For 7 = 1.5, JCii(aj) at this T exhibits 
four sharp modes at w/|Ji| ~ 1,2,3,4, respectively. For 
7 = — 1, there are three broad SCuluj) modes, with rather 
well-defined peaks at w/|Ji| ~ 0,2,3.6 at this T. As for 
the AFM case, in each case 5C22(w) is a single peak at 
one of the non- vanishing SCii{u;) peak positions. 

From an extensive analysis of many lower-T results, 
we found that there are generally four low-T SCii{uj) FM 
mode frequencies, f2i(7), and that (5C22(^) has one non- 
vanishing mode frequency, ^^1(7). These FM pic- 
tured in Fig. m satisfy 

^^i(7)/>^i = |1 - I/7I , for 7 < -1/2, 

= 3, for 7 > -1/2; (32) 

^2(1)/ Ji = 0, for 7 < -1/2, 

= 1 + 27, for 7 > -1/2; (33) 

l]3(7)/^i =2|l-l/7l, for 7<-l/2, 

= |5 - 27I, for 7 > -1/2; (34) 

n4j)/Ji = II-I/7I, for 7< -1/2, 

= 2|1 -7I, for 7 > -1/2. (35) 
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FIG. 13. Plots of the low-T FM mode frequencies Q.i/Ji 
versus 7. The triangle indicates low-T mode scaling, and the 
circle indicates the disappearance of the mode. See text. 

As for the AFM case with 1 7^ 7 > 1/2, for 7 < -1/2, 
the low-T FM mode with frequency ^2 (7) = is a central 
peak, and fJi and U.^ are finite and degenerate. For 7 > 
— 1/2, there are four FM magnon modes, none of which 
is a central peak. The amplitude of the weakest mode 
with frequency ^3(7) is oc T^ as T — > for 7 > —1/2. 
At 7 = 5/2, 1^3(5/2) oc T, so the mode with frequency 
^^3(7) is never a central peak. 

For the low-T FM equilateral triangle, indicated by the 
circle in Fig. |l^, 5C22{^) = ^Cii{ijj), since the amplitude 
of the mode with frequency 1^4(7) vanishes as 7 — > 1, and 
the remaining mode frequencies are degenerate. Hence, 
5Cii{lj) approaches the FM scaling form shown in Fig. ||. 
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FIG. 14. Plots at various low T values of 5Cii(aj) vs. the 
scaled frequency tjlal'^^^/l Ji| for the FM r22(— 0.5) mode. 

We now examine the other special FM case, 7 — —1/2. 
The mode f22(— 1/2), indicated by the triangle in Fig. 
p^ , is a low-T FM mode with scaling properties remark- 
ably similar to those of the AFM mode with frequency 
172(1/2), pictured in Fig. |l^. As T 0, the 5Cii{uj) FM 
mode frequency ri2(— 1/2) \ J2\/\a\^^'^ ^ and the mode 
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shape approaches a scahng function of lal^/^w/l Ji|, as 
pictured in Fig. By comparing Fig. |lj with Fig. |l^, 
the lowest T curves in the two figures are almost identi- 
cal, provided that one rescales T by a factor of 8. Since in 
these special cases, the scaling does not involve scaling of 
the mode amplitude, it does not represent a contribution 
to 5Cii{t) involving a scaling of the time, as for the AFM 
equilateral triangle pictured in Fig. or more generally 
for the AFM low-T equivalent neighbor and four-spin 
ring models. |^,^ It is therefore a new form of mode 
scaling, not previously found in any SMM system. 

V. SUMMARY AND DISCUSSION 

We presented the exact solution for the thermal equi- 
librium dynamics of three classical Heisenberg spins on 
a triangle with two exchange couplings Ji and J2. Al- 
though one might expect that for 7 = J2/J1 ~ 1 the 
behavior would not be too different from that of the 
equilateral triangle, 7 = 1, we found that this is not 
the case, regardless of T. Instead of the two Cij{t) re- 
lated by a sum rule for 7 = 1, for 7 7^ 1 there are four 
Cij{t) related by a sum rule. For 7 0, the piecewise 
continuous Fourier transforms 5Cij{uj) demonstrate the 
profound effects of the absence of translational symme- 
try within the triangle with 7 7^ 1. As T ^ 00, they 
allow us to determine the long-time behavior of the au- 
tocorrelation functions Cii(i) and C22it), which approach 
their distinct long-time asymptotic limits differently, as 
t^^ and t^^, respectively. For the three-spin chain, 7 = 0, 
the autocorrelation function on the chain end, SCu^cu), 
is discontinuous at (jj/\Ji\ = 1 at all T, leading to the 
characteristic t^^ behavior of Cii(i) as T — > 00. 

At low T, setting 7 7^ 1 leads to qualitatively different 
behavior from that obtained for 7 = 1. Regardless of 
the sign of Ji, for 7 7^ 1, there are four low-T modes 
in (5Cii(ijj). For Ji < and 7^1, two of these low- 
T modes, ili and fl^, are degenerate for 7 > 1/2. For 
Ji > 0, and 1^4 are degenerate for 7 < —1/2. In both 
of these FM and AFM regimes, one of the non-degenerate 
low-T modes is a central peak with a width determined 
by the Ji, which grows in intensity as T — > 0. This 
central peak, which can have the largest mode intensity, 
arises from the pinning of a magnon mode in the isosceles 
triangle with Ji 7^ J2, and does not appear in any other 
model which has been solved exactly. There is also a peak 
at Lo — arising from the long-time asymptotic limit of 
Cii{t), which would be broadened by relaxation processes 
not included in our model. 

For the three-spin chain, 7 = 0, the effects of the ab- 
sence of translational symmetry in the chain are even 
more dramatic, leading to qualitatively different dynam- 
ical behaviors of the spins on the chain ends from that of 
the spin at the center. For both signs of Ji, SC22{^) has 
only one mode energy | Ji|. But dCulcu) is distinctly dif- 
ferent. For the AFM three-spin chain, SCii{llj) has three 



low-T modes at n|Ji| with n = 1,2,3, where the n = 1 
mode is doubly degenerate. For the FM three-spin chain, 
SCii{iLj) has four modes at n|Ji| with n — 1,2,3,5. Our 
results suggest that for A^-spin chains, the spins on or 
near the chain ends would have a richer dynamics than 
would those of more central spins, or of spins on closed 
rings of N spins. It would be interesting to see if such 
differences are indeed maintained for larger spin chains 
and for larger spin clusters with two or more Heisenberg 
spin exchange interactions. 

For the four-spin ring with equal near-neighbor ex- 
change couplings Ji, the single autocorrelation function 
6Cii{u!) has a single mode that approaches the fixed fre- 
quency 2|Ji| as T 0, in addition to the AFM (FM) 
mode of the four-spin equivalent neighbor model with a 
frequency that approaches (4| Ji |). (2^j2^ None of these 
is a central peak at T 7^ 0. The isosceles triangle has a 
much richer spectrum of four modes, with tunable fre- 
quencies that depend upon 7, and a central peak for a 
semi-infinite range of 7 values. In addition, at the on- 
set of the central peak, that mode form exhibits low-T 
frequency scaling, instead of the time scaling present in 
the AFM equivalent neighbor and four-spin ring models, 
and it is a scaling function of a differently scaled fre- 
quency than in the FM equivalent neighbor model and 
four-spin ring. (2^,^^ The simplest integrable four-spin 
system with dynamics similar to that of the isosceles tri- 
angle is the squashed tetrahedron, which also involves 
two different near-neighbor exchange couplings. p5| ] 

Inelastic neutron scattering experiments on large sin- 
gle crystals would be a particularly useful technique to 
observe the effects predicted here. By appropriately vary- 
ing the scattering wave vector, all four of the dCij{io) can 
be measured as functions of u; and T. Although experi- 
mental observation of the predicted low-T modes might 
be difficult, since quantum effects are expected to domi- 
nate at very low T, the presence of such modes for three 
classical Heisenberg spins on an isosceles triangle under- 
scores the qualitative changes that occur with different 
exchange couplings. Our numerical results indicate that 
the development of these additional modes appears at T 
values high enough for the classical treatment to be valid. 
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APPENDIX 

We first discuss the T ^ 00 limit of the integrals 
I^{t), Eqs. (|l|)-(|l8|). From Eq. (H), we obtain SC22{t) = 
Ii{t)/A. As T ^ 00, dC22{t) coincides with that of the 
N — 3 equivalent neighbor model, so limx^oo Iijt) is 
obtained by setting A^ = 3 in Eqs. (4)-(9) of Ref. 
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Setting t* = Jit as before, for l2{t) we find, 



lim hit) = / dxf2{x) cos [(1 - 7)0;**] , (36) 

T— >oo 



f2ix) = 



256a;2 



Ax{x^ + 1) - [x^ - l)^ln 



x + 1 
x-1 



(37) 



Then Eq. (||) is derived from /2(0) = /2(2) = 0, and 

Since /2(a;) and f2{x) are continuous at a; = 1, Eq. (^) 
leads to the constants ai and 02 in Eq. ( |23|). The FT of 
/3(t) at arbitrary T is given in Eqs. (|45|)-([48|). 

We illustrate the procedure for computing the FTs by 
obtaining (5C22(w). The FT of cos(sr) is 5{Cj + s) + 5{Cj- 
s), where uj — lo/Ji. For a) > 0, the 6{uj + s) term is 
irrelevant. Then uj satisfies |a; — l|<ci)<a; + l. For fixed 
UJ, the integration interval for x is then determined as a 
fimction of uj. Hence, 

<.min[2,(<i+l)] 

- —6(3 - (i) / dxT22{x), (39) 



AZ 



hai^j) = 6(3 + 2A - C})e{Cj - 1) X 

|.min[2,(Li+l)/(l+A)] 
X / dxTii^3-{x) + 

J(<i-1)/(1+A) 

+6(l-w)6(w- A) X 

^(l+i)/(l+A) 

X / dxTii^3-ix), (45) 

J(l-i)/(l-A) 

/3b(c:') 6(3 - 2A - - 1) X 

X / da;:rii 3+(a;) + 

Jp_l)/(1_A) 

^min[2,(l+ii)/(l-A)] 

+6(1 -ci) / dxTii,3+{x), (46) 

J(l-<i)/(l+A) 

/3c(w) = 6(A - w) X 

<.min[2,(l-tr))/(l-A)] 

X / dx[^Tn.3-{x)], (47) 

J(l+<i)/(l+A) 

(l-xV4)[(cD±a:(l + A))^-l]^ 
(a) ± xA)x^ 

For other values of 7, we obtain similar expressions. 
As r — > cxD, one has 



(48) 



li-il 



lim /2(w) = -r; r 



e(2-cS)/2(cD). 



(49) 



•^22(2;) = e 



_ ^ap" + (7-l)2;"]; 



2(1) 



, (40) 



and 6(z) is Heaviside's step function, 6(z) = if z < 0, 
e(z) 1 if z > 0. 

The exact calculation of SCu is substantially more in- 
volved. We define the functions /a(w) to be the FT's of 
the Ia{t), Eqs. (ra ~ 



The general expression for limT^oo-f3(w) is too com- 
plicated to be given here. For the chain, 

^lim ir'ii^) - 6(5 - ^)eiCu - 3)^4=° (^) 



+eii~cj)l-r,fi^), (50) 



6Cu{uj) = Ii{uj) + hiuj) + /3(^), 
where Ii{uj) — 5C22{'-^)/^- We then find that 



xs 1 - 



1 - 



2sCo 



and UJ — uj/[Ji{l — 7)]. We first write Iz{uj) as 



^3(^) = 77777 ha{uj) + hb{i^) + /3c (tl-) 
16Z L 



For < 7 < 2, setting A = |1 — 7I, we then have 



(41) 



M2-UJ) ds^ii,2(s), (42) 



(43) 



(44) 



p^^o.-x _ (5-w)(-24- 195a;-h229a;2-h47w3_^7a)'') 



^r'^(-) 



A%Cj 

(w^ - 1)(1 -f 7(I'2) CZ--1 



In- 



[Co^ - a){Cj^ - ly ^ c^ + i 



In 



4J)2 

iCo^ + 67^2 _ 24 

^"^^ e;^ 

',-,2 iU,~,4 I OQ,~,2 



2(c:; - 2)' 



(d)^ - l)(t:'4 + 23c;;^ + 8) tl- + 1 



4(1-2 

7(1;'* - 66^2 + 951 
^3c K^) = 



' UJ-l' 



1 -(2)2 



(51) 



(52) 



{1-Cj^){1 + 7Cj^) ^ l-uj'' 

uj'^ 16 
(4-.2)(,_,2)2^^4(4_^^ (53) 
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The leading asymptotic behavior, Eq. (E2|), arises from [9 



4r°(5) 


= 0, 




(54) 


[10] 


/3r°(3) 




592, „ 
--ln2, 


(55) 


[11] 




= 4^(1)- 


119 


(56) 


[12] 


4r°(o) 


^H^_331n2. 

n 




(57) 


[13] 



We now derive the AFM low-T mode frequency ^^2(7)- 
As a — > — 00, the integrand in Eq. ( |39| ) becomes sharply 
peaked about some w-dependent value xq. We then per- 
form an asymptotic evaluation of the form, 



J a 



(Ph 



(58) 



Xq 



h{x) = a [o;^ + (7 - \)x^] + Incj + 
(l+c:;2„^2)2 



In 



1 - 



4a)' 



r,2 



(59) 



Inside the integration region, the integrand has exactly 
one maximum at xo- For 7 > 1, x^^''^^ = |1 — cij|+e, while 

for 7 < 1, xj,''^^^ = 1+tJ-e, where e > is 0{a^'^). We 
then evaluate the integral in Eq. (|5^). It is maximal for 
Lo — where 



f^i(7)/|Ji| = |l-l/7l 



(60) 



For 7 < 1/2, this result gives a spurious maximum, since 
it would require x > 2. Since the result of the asymptotic 
evaluation is a monotonically increasing function of lo for 
< 1, we therefore conclude that fii — \ Ji \ for 7 < 1/2. 
The proof of Eq. (pq) is thus complete. 
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